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Abstract. We adapt and modify the eigenfunction method of computing the power- 
law spectrum of particles accelerated at a relativistic shock front via the first-order 
• Fermi process [6] to apply to shocks of arbitrarily high Lorentz factor. The power-law 

| index of accelerated particles undergoing isotropic small-angle scattering at an ultra- 

relativistic, unmagnetized shock is found to be s = 4.23 ± 0.2 (where s = din f /d\np, 
with / the Lorentz- invariant phase-space density and p the momentum), in agreement 
with the results of Monte-Carlo simulations. We present results for shocks in plasmas 
with different equations of state and for Lorentz factors ranging from 5 to infinity. 
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We study a stationary shock front in the x — y— plane. The accelerated particles 
are assumed to be test-particles without influence on the dynamics of the plasma 
or the jump conditions at the shock-front. The plasma flows along the 2-axis, with 
constant velocities w_ in the upstream (z < 0) region and u + downstream (z > 0), 
the velocities are related by the Rankine-Hugoniot jump conditions. 

Test-particles are injected into the acceleration process and their interaction with 
the plasma flow is assumed to give rise to diffusion in the angle cos" 1 fi between a 
particle's velocity and the shock normal. In the frame of the shock front this leads 
to a stationary transport equation valid for the local plasma rest frame and given 
in mixed coordinates as [6] 



r(« + /% = ^(i-M^ (i) 



where the plasma speed u is measured in units of the speed of light, T = (1 — u 2 ) 
is the Lorentz-factor, f(p,fi,z) is the (Lorentz invariant) phase-space density as 



homepage: http://www.mpi-hd.mpg.de/theory/ 



a function of the particle momentum p, direction p and position, p and p are 
measured in the local rest frame of the plasma, whereas z is measured in the rest 
frame of the shock front. 

Equation (1) is solved using the separation Ansatz [6] 

+00 

f(p,u,p,x) = 9i(p)Qi(n,u)exp(Aiz/r) , (2) 

i=— 00 

valid in each half-plane with Aj and Qi the eigenvalues and eigenfunctions of the 
equation 

-A i (u + fi)^Q i ((i,u) = (3) 

The momentum distribution of particles with energy far above the injection energy 
range - those in which we are interested - takes the shape of a power-law gi(p) oc p~ s 
with a power-law index s, since there is no preferred momentum scale in this range. 

Matching the expansion (2) across the shock front according to Liouville's The- 
orem and imposing physically realistic boundary conditions up and downstream 
leads to a nonlinear algebraic equation for the power law index s. 

In [6] and [4] only the eigenfunctions with % < were used and the method was 
applied to mildly relativistic shock speeds (r_ < 5). Here, we use the eigenfunctions 
with % > and calculate them directly with a numerical scheme. In the limit w_ — > 1 
an analytic expression is available [7]. Four eigenfunctions (i = 1,3,5, 7) are shown 
in Fig. 1A as functions of the cosine p s = (p + u)/(l + pu) of the angle between 
the particle direction and the shock normal, measured in the shock rest frame. For 
% > 1 they are oscillatory for — 1 < p s < and for all i > fall off monotonically 
in the range < p s < 1. 
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RESULTS 

The index s of the momentum spectra of the accelerated particles in different 
cases are shown in Fig. IB. The jump conditions investigated are those for a rel- 
ativistic gas both up and downstream: u_u + = 1/3 and for a strong shock in a 
medium with adiabatic index 4/3 [5]. 

Also we investigate two different scattering operators, D^i = 1 — p 2 (isotropic 
small/angle scattering) and = (1 — p 2 ) x (p 2 + O.OI) 1 / 3 corresponding to 

scattering in weak Kolmogorov turbulence, together with a rough prescription for 
avoiding the lack of scattering at p — [4] . For high upstream Lorentz-factors the 
power-law index settles at a value around 4.23 for all equations of state, which is 
reproduced in the limiting case u_ — > 1. The scattering operator has only a minor 
effect. 
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FIGURE 1. A) (left) The eigenfunctions Q, for i = 1, 3, 5, 7 for T_ = 223, as a function of the 
(cosine of the) angle between the particle speed and the shock normal, measured in the shock 
frame, for a relativistic gas. B) (right) The power-law index s for relativistic gas with isotropic 
(solid, 2nd from top) and anisotropic (dashed-dotted, top) scattering operator and for a strong 
shock in a gas of adiabatic index 4/3, with isotropic (dashed, 4th from top) and anisotropic 
(dotted, 3rd from top) scattering 



These results are in agreement with the asymptotic Monte-Carlo results of Gal- 
lant et al. [3] and those of Bednarz & Ostrowski [1] for T_ ps 200. Anisotropic 
scattering, which has not been treated by Monte-Carlo simulations, leads to a slight 
steeping in the power-law spectrum, because fewer particles are able to cross the re- 
gion ji ps and return to the shock. From observations of GRB afterglows, Galama 
et al. [2] and Waxman [8] have found synchrotron spectral indices corresponding 
to s ~ 4.25, implying that the particles could indeed have been accelerated by the 
first order Fermi mechanism operating at an ultrarelativistic shock front. 
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